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Abstract 

Gravity can be thought as an emergent phenomenon and it has a nice "thermodynamic" 
structure. In this context, it is then possible to study the thermodynamics without knowing 
the details of the underlying microscopic degrees of freedom. Here, based on the ordinary 
thermodynamics, we investigate the phase transition of the static, spherically symmetric 
charged black hole solution with arbitrary scalar curvature 2k in Horava-Lifshitz gravity at 
the Lifshitz point z = 3. The analysis is done using the canonical ensemble frame work; i.e. 
the charge is kept fixed. We find (a) for both k — and k = 1, there is no phase transition, 
(b) while k — — 1 case exhibits the second order phase transition within the physical region 
of the black hole. The critical point of second order phase transition is obtained by the 
divergence of the heat capacity at constant charge. Near the critical point, we find the 
various critical exponents. It is also observed that they sati sfy the usual thermodynamic 
scaling laws. 



1 Introduction 

One of the fascinating solutions of the gravity theories is the black hole space time. Pioneering 
works of Bekenstain [T], Hawking [2] and later on Bardeen et al [3] illuminated the fact that 
all the laws of black hole mechanics are identical to those of the ordinary thermodynamics with 
properly identification of the respective quantities, like temperature, entropy, energy etc. and 
hence they do act as the ordinary thermal system. People then started thinking that black holes 
may play a crucial role to give idea on the quantum nature of the gravity. But till now there is 
no such complete theory of quantum gravity. 

Several recent results strongly indicate the possibility that the field equations of gravity have 
the same status as the equations of fluid mechanics or elasticity. For a recent review, see [I]. One 
specific implementation of this idea considers the field equations of the theory to be 'emergent' 
in a well-defined sense, rather than use that term in a more speculative vein like e.g., considering 
the space and time themselves to be emergent etc. The evidence for such a specific interpretation 
comes from different facts like the possibility of interpreting the field equation in a wide class of 
theories as thermodynamic relations [5] , the nature of action functional in gravitational theories 
and their thermodynamic interpretation [B] , the possibility of obtaining the field equations from 
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a thermodynamic extremum principle [7] , application of equipartition ideas to obtain the density 
of microscopic degrees of freedom [8] , the equivalence of Einstein's field equations to the Navier- 
Stokes equations near a null surface [9] etc. The important fact in the emergent paradigm is 
that one does not need to know about the details of the microscopic description of the theory. 
These lead to the idea that one can study the different aspects of thermodynamics, in the case 
of gravity, without any detail of the underlying microscopic structure. Here we will adopt the 
same logic and study the phase transition of the black hole solution in Hof ava-Lifshitz theory 
following the prescription in ordinary thermodynamics. 

Studying the thermodynamics |10|-|27| as well as the critical behavior |28]-|41j of various 
black hole solutions in the framework of usual Einstein gravity has been a fascinating topic of 
research for the past few decades. Apart from these, there also exist topological black hole 
solutions whose thermodynamics has been studied in Einstein gravity \42\ I43j. Einstein- Gauss- 
Bonnet gravity and dilaton gravity |44j as well as Lovelock gravity [45 i. On the other hand, 
inspired by the dynamical critical phenomena in usual condensed matter systems, very recently 
P. Hofava proposed a UV complete theory of gravity [36] that reduces to the usual Einstein 
gravity at large scales. Since then a number of attempts have been made in order to understand 
various aspects of this theory [37]-[lH] including different cosmological aspects [50J. Various 
black hole solutions were also found in [JH E21 E3] • The thermodynamics of these black holes 
have been studied in [52j [5H [55]. Although these attempts are self contained and rigorous, 
still there remains some major questions which have not yet been attempted. Studying the 
critical behavior of black holes is one of the important aspects, which we aim to explore for 
Hofava-Lifshitz theory of gravity. 

In ordinary thermodynamics, critical exponents plays a crucial role in order to understand 
the singular behavior of various thermodynamic entities near the critical point. In usual ther- 
modynamic systems, it is customary to express the singular behavior of various thermodynamic 
systems in terms of power laws characterized by a set of static critical exponents which deter- 
mines the qualitative nature of the phase transition near the critical point. All these exponents 
are not independent and are found to satisfy certain thermodynamic scaling laws near the crit- 
ical point |56]-|57j. It is generally observed that for d > 4, where d is the spatial dimension of 
the system, the critical exponents do not depend on the spatial dimension of the system, which 
reflects the mean field characteristics of a theory. On the other hand if the interaction between 
the constituent elements are short range then the critical exponents are found to be depe ndent 
on d. 

In this paper, based on a canonical framework (i.e; keeping the charge (Q) of the black hole 
fixed [13J), we investigate the critical behavior of topological charged black holes in Hofava- 
Lifshitz theory of gravity at the Lifshitz point z = 3. The black hole solution is given in |52| . 
Although all the thermodynamic quantities were evaluated earlier [52], a detailed study of the 
nature of phase transition is still lacking. Particularly the issue regarding the scaling behavior 
of (charged) black holes has never been investigated so far in the frame work of Hofava-Lifshitz 
theory of gravity. We investigate the phase transition phenomena for all the three cases taking 
k = 0, ±1. We observe the following interesting features: 

• There is no Hawking Page transition [10] for black hole with k = 0, 1. 

• For k = —1, there is a upper bound in the value of the event horizon, above which the 
temperature becomes negative. This indicates that above this critical value of horizon, the 
black hole solution does not exist. We will call this valid range as the physical region. 

• Within the physical region, interesting phase structure could be observed for the hyperbolic 
charged black holes (k = —1). In this particular case we observe the second order transition. 

This is different from the usual one. In Einstein gravity, Hawking page transition occurs for k = 1 
while there is no such transition for k = 0,-1. The critical point of second order phase transition 
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is marked by the divergence in the heat capacity at constant charge (Cq). Finally, we explicitly 
calculate all the static critical exponents associated with the second order transition, and check 
the validity of thermodynamic scaling laws near the critical point. Interestingly enough it is 
found that the hyperbolic charged black holes (k = — 1) in the Hofava-Lifshitz theory of gravity 
fall under the same universality class to that with the black holes having spherically symmetric 
topology (k = 1) in the usual Einstein gravity. The values of these critical exponents indeed 
suggest a universal mean field behavior in black holes which valid in both the Einstein as well 
as Hofava-Lifshitz theory of gravity [40]-|41j. 

The plan of the paper is as follows: We begin in section 2 by giving a brief introduction of 
the black hole solution in Hofava-Lifshitz gravity and use it in section 3 to study the different 
thermodynamic quantities as well as the phase transition. In section 4, the critical exponents 
near the critical point are being evaluated and a brief discussion on the validity of the ordinary 
scaling laws is presented. Finally, we conclude in section 5. 



2 Charged black hole space-time in Hofava-Lifshitz gravity 

In this section a brief discussion on the black hole solution in Hofava-Lifshitz gravity will be 
presented. For details, one can follow [52] where the meaning of all the parameters are given 
explicitly. We will mainly concentrate on the solutions at the Lifshitz point z = 3, particularly, 
the static spherically symmetric topological charged black holes solution 

dr 2 

d S 2 = -f(r)dt 2 + — + r 2 dfl 2 , (1) 

where, dQ 2 , is the line element for a two dimensional Einstein space with constant scalar curvature 
2k. Without loss of generality, one can take k = 0, ±1 respectively. The form of the metric 
coefficient, for the detailed balance condition, is given by [52] 



f(r) = k + x 2 - \jux-^ , (2) 

where x = yf—Kr and A(= — Jx) corresponds to the negative cosmological constant. The physical 
mass and the charge (Q) corresponding to the black hole solution are respectively given by, 

K 2 fi 2 n k V^X K 2 fi 2 n k V^A 

M = a: Q = q (3) 

16 16 K ' 

where, a, q are the integration constants, is the volume of the two dimensional Einstein 
space and k, fi are the constant parameters of the theory. The event horizon is the solution of 
the equation f(r+) = 0. 



3 Phase transition 



In ordinary thermodynamics the phase transition is studied by the divergence of relevant ther- 
modynamic quantities. Here, the same technique will be adopted. We shall first calculate the 
Hawking temperature, entropy and specific heat of the black hole using a canonical ensemble 
frame work, which means that we shall carry out our analysis keeping the total charge (Q) of 
the black hole fixed |13j . Finally, a graphical analysis will be given to study the phase transition. 
Before we proceed further, let us mention that in the following analysis we re-scale our variables 
as M -> S -)• #- , Q -> Q / - x , where we have set 16vrA = 1 with A = 
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The Hawking temperature is calculated as 



_ f'(r + ) _ V^A(3x 4 + + 2kx\ - fc 2 - f) 

4tt 8irx + (k + x 2 + ) ' 1 ' 

and after the above mentioned re-scaling the entropy is found to be 

_ fdM (x\ k \ 

S = J — = ( -f + -lnx + \ + So . (5) 

In the above, So is the integration constant which must be fixed by the physical consideration. 
To be precise, entropy is always determined upto some additive constant. But in all physical 
considerations the difference is important. Therefore in the present analysis we will consider 
only the difference. However, the integration constant So could be determined using the usual 
thermodynamic prescription; i.e. determination of entropy in the T — > limit, which might be 
the entropy of an extremal black hole. 




Figure 1: Temperature (T) plot of topological black holes for k = 0, 1 with respect to x+ for 
Q = 5 and 1 = 1. 




Figure 2: Temperature (T) and Entropy (S) plot of topological black holes for k = — 1 with 
respect to x+ for Q = 5 and 1 = 1. 

Before we proceed further, let us first try to analyze the behavior of Hawking temperature 
(T) for different choices of k. First of all, for the case k = 0, 1, we note that the Hawking 
temperature (T) is a monotonically increasing function of horizon radius (see fig. 1) which 
indicates that the corresponding heat capacities (Cq) are always positive definite. Therefore 
these black holes are globally stable and there is no phase transition. 

Next consider the other case, k = —1. For this, the Hawking temperature Q reduces to the 
following form: 



It it interesting to note that the Hawking temperature (T) is always positive in the range 
< x+ < 1. Furthermore, if we plot the entropy (S) as a function of x+, using ([5]), it shows 
that S is also positive in this range (see fig. 2). Interestingly enough this condition is found 
to be valid for any value of the physical charge (Q) of the black hole. On the other hand, we 
note that T — > — oo as x+ —> 1. This is due to the fact that as x+ — > 1 the second term on the 
r.h.s of ^ dominates over the first one which ultimately produces a large negative temperature. 
This is a nonphysical situation and the corresponding black hole solution does not exist for 
x + > 1. Considering this fact, in the present paper we carry out our analysis in the physical 
range < x+ < 1 where the temperature (T) of the black hole is finite as well as positive 
definite. 

In the above specified range (0 < x+ < 1) we observe a change in slope at x+ = xi of the 
corresponding (T — x + ) plot (see fig. 2). This change in slope signals a discontinuity at x + = X2 
in the corresponding heat capacity Cq. It gives a indication of second order phase transition 
at x+ = X2, which we shall refer as the critical point of phase transition. In order to make the 
discussion more transparent we next compute the corresponding heat capacity (Cq). 

Using Q and ([5]) the heat capacity is determined as 

( M. 



(k + 4) 2 (34 + 2k.vi - k 2 - ( 4 ) 



64 + 14fc4 + 10/c 2 4 + 2 A; 3 + Q 2 (k + 3x 2 




Figure 3: Specific heat (Cq) plot of topological black holes for k = —1 with respect to x+ for 
Q = 5 and 1 = 1. 

From the above figure (fig. 3), we see that specific heat (Cq) indeed suffers a discontinuity at 
x + = x%. This shows that there is a genuine second order phase transition at x+ = xi- For the 
phase 1 the heat capacity (Cq) is always found to be positive (see fig. 3) which means that this 
phase is thermodynamically stable. On the other hand, Cq < for the phase 2 and therefore it 
is an unstable phase. 

This scenario is completely different in the Einstein gravity, where we have no Hawking Page 
transition for k = 0,-1. Whereas, k = 1 case exhibits Hawking Page transition |42j-[44j. For 
the above phase structure one can further investigate the critical behavior of the black hole in 
Hofava-Lifshitz theory of gravity near the critical point X2, which we aim to discuss in the next 
section. 
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4 Critical phenomena and scaling laws 



In this section, based on a thermodynamic approach, we explicitly calculate various critical 
exponents (that is associated with the second order phase transition at x 2 ) for the hyperbolic 
charged black holes (with k = — 1) in the fixed charge (Q) ensemble. 

In the theory of phase transitions, critical exponents play an important role in order to 
understand the singular behavior of various thermodynamic entities near the critical point(s). 
In order to have a complete understanding of the physics of phase transition phenomena, one 
generally introduces a set of static critical exponents (a, (3, 7, 5, if, 1/1, v, rf) which play a central 
role in the theory of critical phenomena, and studies the so called thermodynamic scaling laws 
near the critical point. 

In order to find the critical exponents, we write near the critical point 



x + = x 2 (l + A), 



(8) 



where |A| << 1. Defining T 2 = T(x 2 ), the Taylor expansion of the temperature (T) about 
x + = X2 yields 



T 2 + 



dT 

dx+ 



Q 



-x 2 ) + - 



1 



X + =X 2 



d 2 T 

dx\ 



X + =X 2 



{x + —x 2 ) 2 + higher order terms . (9) 



It has been shown earlier that at the critical point x 2 , Cq diverges (see Figure 3). Therefore 
implies that the second term on the right hand side of ^ vanishes. Hence neglecting the 
higher order terms and then using Q in ^ we obtain, 



A 



2 (T-T 2 ) 

X2 



1/2 



where, 



D 



d 2 T 

OX 2 : 



6x'i + 3Q 2 xj - <ox\ - <oQ l x\ -2-Q 



2^4 
2 



4(4 - 1) 3 



X + =X 2 

Now to find the critical exponent a which is defined by the standard relation 

~ |i — J-2\ , 

we use the relation Q in ([7]) and then keep the terms only linear in A to obtain 

M(x 2 ,Q) 



A(36x^ - 56x| + 20x| + 6Q 2 x^) ' 



with 



N{x 2 , Q) = {x\ - l) 2 (3x| - 2x 2 2 ~ 1 - Q 2 /2) 



Finally, using (10) in the above we find the behavior of Cq near the critical point x+ = x 2 - 

A(x 2 ,Q) 



Cq 



where, 



A(x 2 ,Q) 



(T-T 2 )V2 
N(x 2 ,Q)VDx 2 



V2(36x| - 56x| + 20xj + ®Q 2 4) 



(10) 
(11) 
(12) 

(13) 
(14) 

(15) 

(16) 



A comparison of (15) with the standard relation (12) yields a = 1/2. 
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Next the determination of the critical exponent /3 which is defined through the relation (for 
a fixed value of charge Q) as, 

$(x+)-$(x 2 )~ \T-T 2 f (17) 
will be done. Here the potential $(x) is given by 

Q 



<£(x) 

To proceed, let us first Taylor expand 3>(x + ) = Q- about x + = x 2 : 



$(x + ) = $(x 2 ) + 



dx+J Q 



X+=X2 



(x+ — x 2 ) + higher order terms 



#(x 2 ) — "^(^+ — £2) + higher order terms . 



(18) 



(19) 



As before, neglecting the higher order terms and then using (Tsb we obtain, 



^x + )-^{x 2 ) = -%M(T-T 2 f/\ 



(20) 



where in the final step (10) also has been used. This immediately determines (3 = 1/2. 

To find out the exponent 7, associated with the divergence of the inverse of the isothermal 



compressibility K T = Q 



8Q 



32], we first find the explicit expression of K~ x . The 



thermodynamic identity (gf)g(^j (gi 



T yields 



VdQJT \dT)Q\dQ)<s>' 



/d<S>\ /dT\ 



(21) 



In order to evaluate the right hand side of (21) first note that, 



Also, from the functional relation, 
we find, 



8<s> 

\df/Q far_ 

8x+ 



T = T(x + ,Q) 



(22) 



(23) 
(24) 



/dT\ ( dT \ ( dx +\ ( dT \ 
VdQJ* ~ \8x^)q ' '* + \dQJ*+' 

Right hand side of both ( 22 ) and ( 24 ) can be easily calculated using the relations Q and ( |18[ ) 
which finally yields, 



K T X = Q 



d<P\ 
dQJ, 



Q \ 6x° + Wx% - Ux% - 2 + Q\x% + 1) 
x~l) §x\ + 10x2 - 14x1 - 2 + Q 2 (3x2_ - i; 



(25) 



Then to obtain the near critical point expression, use (|8j) and next (10) in the above. This 
yields, 

R -i „ QVDjQxl + lOxj - 14x| - 2 + Q 2 (x| + 1)) _ 1/2 

T v/2 (36x| _ 5 g x 4 + 20x 2 + 6 Q2 X 2) 2 1 J 

Comparing with the standard relation ~ |T — T 2 |~ 7 , defined for the fixed value of charge 
(Q), we get 7 = 1/2. 
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The critical exponent 5 is denned through the relation, 



ii/i 



(27) 



at constant temperature T. To find it, we first expand Q(x+) in a sufficiently small neighborhood 
of = X2 which yields, 



\dx + J T 



(x + - x 2 ) + 



1 



x+=a;2 



Q(x+) = Q(x 2 ) + 

+ higher order terms . 
Since T = T(x+,Q) (see Eq. @), one finds for the fixed T 



\dx\J T 



(x + - x 2 ) 



x + =x 2 



(28) 



dQ 
dx+ 



dT \ (dQ 
dx+ J n \ dT 



Now at the critical point x+ = x 2 , Cq diverges and so, as earlier, 



dT 

dx+ 



Q 



x+=x 2 



fore, after neglecting the higher order terms, (28) reduces to 

T 



x + - x 2 



{Q{x + )-Q{x 2 )) 1 ' 2 



where, 



M 



dx\ 



x + =x 2 



18x1 + S2x 2 ~ Ux 2 ~ 3Q 2 x% -Ax 2 2 -2-Q 2 
2Q 2 x 2 (x 2 - l) 2 ' 



(29) 
0. There- 

(30) 

(31) 



Next since $>(x+) = J*-, we have the functional relation $ = Q) and so 



d<S> 



dx+ ) t \dx 



<9<I» 



+ 



d$\ f dQ 



dQ/x+ \dx + J t 



Therefore using (29) we find, 



iff ) ' 




[(f ) 




\dx+J T 


X+=X2 


\dx + J 





Q 



,-2 ' 



(32) 



(33) 



X + =X2 



Finally, expanding close to the critical point x + ~ x 2 at constant T and then using (30) 

and (33) we obtain, 



Q(sc + ) - Q(x 2 ] 



1/2 



(34) 



Now re-expressing (15) by using (10) in the following form: 

Y A 12 A 



Hence, the critical exponent is read off as 5 = 2 

■ usi 

Cq ~ V D x 2 A ~ V D (x+ - x 2 ) 1 
and then substituting the value of x + — x 2 from ( 30 ) we obtain^ 

1 



Cq 



(Q( X+ ) - g 2 )V2 • 



(35) 



(36) 



1 Q(x 2 ) = Qa- 
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Comparing this with the standard definition Cq ~ \Q(x + ) — Q 2 \ v we find ip = 1/2. 
In the following the calculation, the critical exponent ip, defined by 

S(x+)-S(x 2 ) ~ \Q-Q 2 f , (37) 

will be found out. Expansion of the entropy for fixed charge Q about the critical point x 2 yields 

S(x + ) = S(x 2 ) + ( S—) (%+ — x 2 ) + higher order terms . (38) 

L\ax+/ Qix + =x 2 



Then following the identical steps as earlier and using ( 30 ) we find 



8(x+) - S(x 2 ) ~ M-Hyi(Q - Q 2 )V2 ; (39) 
which yields ip = 1/2. 

In the following table, we give the values of the critical exponents for the present example: It 
Table 1: Various critical exponents and their values 



Critical Exponents 


a f3 


7 5 


f 




Values 


1/2 1/2 


1/2 2 


1/2 


1/2 



may be interesting to mention that the critical exponents for the second order phase transition 
in the case of hyperbolic black hole in Hofava-Lifshitz gravity are exactly equal to those obtained 
earlier in [3D] for the Einstein-Born- Infeld gravity with k = 1. 

In the last part of the section, we discuss about the validity of the usual thermodynamic 
scaling laws for our present situation. It should be mentioned that in the realistic thermody- 
namic systems various critical exponents are found to satisfy certain (scaling) relations among 
themselves, known as thermodynamic scaling laws which may be expressed as, 

a + 2/3 + 7 = 2, a + (3(5 + 1) = 2, (2 - a)(6ip - 1) + 1 = (1 - a)6 

j(6+l) = (2-a)(6-l), 7 = /3(S-l), ^ + 2^-6~ 1 = 1. (40) 

Since the critical exponents, we find here, are exactly equal to those obtained in the usual 



Einstein-Born- Infeld gravity for k = 1 [10], all the above relations (40) are indeed satisfied for 



the present case, in spite of difference in the nature of the black hole solution and the phase 



structure. This can be checked easily by substituting the values from the table in (40). In this 
sense, it may be possible that the black holes in both type of theories fall under the same class 
such that their critical exponents are identical. 

Finally we are in a position to check the additional scaling laws for the exponents v and n which 
are associated with the diverging nature of correlation length and correlation function near the 
critical point. Since the spatial dimension (d) of our theory is < 4, therefore it will be natural 
to assume that the additional scaling relations, 

7 = 1/(2-7/), 2-a = vd (41) 

will hold in general. This immediately determines the other exponents. Taking d = 3 and using 
the exponent values from table 1, we finally obtain, 

u = l/2, 77 = 1. (42) 



In practical situation, the assumption that the scaling laws (41) are valid, may not be true. 
Hence it must be checked that these are indeed valid or not. Alternatively, one should find the 
above exponents by some alternative method. 
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5 Conclusions 



It is now evident that gravity and thermodynamics are closely connected to each other [H El [3] . 
The repeated failure to quantize gravity led to a parallel development where gravity is believed 
to be an emergent phenomenon just like thermodynamics and hydrodynamics instead treating it 
as a fundamental force |4J. The fundamental role of gravity is replaced by thermodynamic inter- 
pretations leading to similar or equivalent results without knowing the underlying microscopic 
details. 

In this paper, we adopted the standard thermodynamic approach to explore the phase struc- 
ture of the topological charged black holes in Hofava-Lifshitz gravity for the Lifshitz point z = 3. 
These black holes were introduced earlier in [52] . In [52] the authors have computed the general 
expression for the Hawking temperature for these black holes. In spite of this attempt, till date 
the issue regarding the nature of phase transition, particularly the issue of critical phenomena, 
remains completely unexplored. In the present paper, we attempted to provide an answer to all 
these questions. A number of interesting features have been observed in this regard which have 
never been explored so far to the best of our knowledge. These are as follows: 

• We showed that there exits a physical range < x + < 1, where the black hole solution exists. 

• Within this range, the second order phase transition occurs. 

• The situation is not exactly identical to Einstein theory. In Einstein theory, there is no 
Hawking Page phase transition for k = 0,-1 while k = 1 case exhibits such phase transition. 
Here, only k = — 1 solution has phase transition. 

• Finally, the critical exponents near the critical point of the phase transition were derived. This 
was done again using the ordinary thermodynamic analogy. It may be noted that the critical 
point has been marked by the discontinuity of the heat capacity (Cq) which suggests that it is 
a second order phase transition. Moreover, as a point of conformation, we also found that the 
critical exponent, associated with the divergence in Cq, is 1/2 which simultaneously indicates 
a mean field feature as well as a second order nature of the phase transition. Interestingly, 
these are exactly equal to those of Einstein-Born- Infeld theory |40]-|41] and also they satisfy the 
ordinary thermodynamic scaling laws. This in fact suggests two remarkable facts: (1) black holes 
in both the Einstein as well as Hofava-Lifshitz theory of gravity fall under the same universality 
class, (2) there exists a universal mean field behavior in both of these gravi ty theories. 

It may be pointed out that the issues, dealt with the present paper, are notoriously difficult 
to study. So far what we have done is completely based on the ordinary thermodynamic analogy. 
But a simple analogy with the thermodynamics may not be enough to draw safe conclusions. 
Also the understanding of the thermodynamic interpretation of gravity is far from complete since 
the arguments are more heuristic than concrete and depend upon specific ansatz or assumptions. 
Therefore the results are suggestive rather than definitive. Hence it is necessary to compute the 
critical exponents by other alternative/direct procedures; e.g. correlation function technique 
of the scalar modes, the AdS/CFT correspondence, re-normalization group scheme, etc. which 
is beyond the scope of the present paper. In spite of this, we still believe that the analysis, 
presented here, could illuminate the underlying microscopic structure of the gravity. 
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